NORMAL AND ANOMALOUS SCALING OF THE FOURTH-ORDER 
CORRELATION FUNCTION OF A RANDOMLY ADVECTED PASSIVE SCALAR 
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For a delta-correlated velocity field, simultaneous correlation functions of a passive scalar sat- 
isfy closed equations. We analyze the equation for the four-point function. To describe a solution 
completely, one has to solve the matching problems at the scale of the source and at the diffusion 
scale. We solve both the matching problems and thus find the dependence of the four-point correla- 
tion function on the diffusion and pumping scale for large space dimensionality d. It is shown that 
anomalous scaling appears in the first order of 1/d perturbation theory. Anomalous dimensions are 
found analytically both for the scalar field and for it's derivatives, in particular, for the dissipation 
field. 
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INTRODUCTION 



It is already a commonplace to talk about an anoma- 
lous scaling of the high-order correlation functions in de- 
veloped turbulence. By this, the deviation of the scaling 
exponents from their "naive" values taken from dimen- 
sional estimate or perturbation theory is usually meant. 
Another meaning ascribed to that term is related to the 
cases where the exponent of the 2n-th correlation func- 
tion is not n times the exponent of the second one so 
that the degree of non-Gaussianity depends on scale. 
The experimental evidence in favor of anomalous scaling 
of a scalar field advected by turbulence exists for some 
time [0-|4| while the attempts of the consistent deriva- 
tion (starting from the equations of fluid mechanics) of 
the correlation functions of the order higher than two 
started only recently |5|-|l0|l. The problem of a passive 
scalar advection, being of physical importance by itself, 
may serve also as a starting point in studying anomalous 
scaling in turbulence At a first step, we found 

||i"2f the whole set of the simultaneous correlation func- 
tions for the Batchelor-Kraichnan problem of a scalar 
advected by a large-scale random velocity field 
It has been shown that, whatever be the (finite) tempo- 
ral correlations of the velocity field, all correlation func- 
tions of the scalar are integer powers of a logarithm for 
all the distances in the convective interval of scales and 
no anomalous scaling thus appears at the leading terms. 
The present paper is an account of the next step: we 
consider multi-scale velocity field with power spectrum. 
Following Kraichnan p5[ , we restrict ourselves by the 
simplest possible temporal behavior assuming both ve- 
locity field and scalar source to be white in time. That 
leads to a substantial simplification of the analytical de- 
scription since any simultaneous correlation function of 
a scalar satisfies closed linear differential equation of the 



second order (see Q and below). 

In an isotropic turbulence, n-point correlation function 
depends on n(n — l)/2 distances for a dimensionality of 
space d > n — 2. For the pair correlation function, the 
respective ordinary differential equation could be readily 
solved for any distance between the points — see jl5| and 



(1.19) below. The solution is expressed via the flux of a 
squared scalar P 2 and molecular and eddy diffusivity, the 
scaling exponent £2 in the inertial interval being fixed by 
the condition of the flux constancy. If the scaling expo- 
nent of the 2n-th correlation function is nC,2 , that is called 
a normal scaling. An anomalous scaling would mean that 
the true answer has additional factors like {l/r) A where 
A is anomalous dimension, r is the distance between 
points and I is some length parameter. One may imag- 
ine three reasons for anomalous scaling according to the 
three lengths that may be relevant: i) I is the pumping 
scale L so that the anomalous scaling appears due to in- 
frared non-locality, ii) I is the diffusion scale so that the 
anomalous scaling appears due to ultraviolet non-locality 
|U|[uJ and iii) anomalous scaling appears due to the exis- 
tence of the high-order integrals of motion so that the 
length parameter appears from the ratio of the different 
fluxes. First nontrivial object that may reveal anoma- 
lous scaling is the fourth-order correlation function. One 
may be interested in two-point objects like ((#1 — (?2) 4 ) 
or pair correlation function (e^) of the dissipation field 
e(t, r) = K[V6*(i,r)] 2 yet they do not satisfy any closed 
equation. To find such two-point-fourth-order objects, 
one should solve the complete equation for the four-point 
function and then fuse some points. Considering d > 2, 
one has to deal with the space of six variables which 
makes the direct solution at arbitrary parameters quite 
difficult. Gawedski and Kupiainen recently devel- 
oped the perturbation theory that describes the effect of 
the weak advection on the fourth-order correlation func- 
tion for a non-cascade diffusion-like regime and found the 
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respective anomalous exponents. Our main target here is 
the description of the cascade in the convective interval 
of scales, we solve the equation for the four-point corre- 
lation function assuming the space dimensionality d to 
be large. Non-gaussian part of the correlation function 
is small in parameter 1/d which makes it possible to de- 
velop a regular perturbation theory. In the present paper, 
we obtain the scaling exponents analytically in the first 
order in 1/d (assuming 1/d to be the smallest parameter, 
in particular, 1/d <C C2) and show that the main contri- 
bution into the general fourth-order correlation function 
has L-related anomalous scaling with 



A = 



4(2 - C2) 



(0.1) 



Note that A turns into zero when (2 — * 2 which con- 
tradicts Kraichnan's closure HQ and is in a qualitative 
agreement with Gawedski and Kupiainen result fl0|| even 
though they considered a different limit. 

For the correlation functions between the points sep- 
arated by some distances from the convective interval 
and by others from a diffusion interval of scales, we show 
that the main term does not depend on r^- We also 
show that the sublcading terms in the case of some dis- 
tances being much smaller than others necessarily have 
another anomalous exponents with respect to small dis- 
tances. Such subleading terms determine, in particu- 
lar, correlations of the spatial derivatives of the scalar 
field, which thus have an r^-related anomalous scaling 
(in addition to a general L-related scaling). Note that 
the results on r^-dependence are obtained without 1/d- 
expansion so they are exact at any dimension. For ex- 
ample, the dissipation field has an r^-related anomalous 
scaling with the exponent that is equal to the scaling 
exponent of the eddy diffusivity so that the irreducible 
correlator ((eie 2 )) oc (L/Yi 2 ) A r° does not depend on the 
diffusion scale when the distance ri 2 is in the con- 
vective interval. Another consequence of our results is 
the statement that (e 2 )/ (e) 2 tends to infinity with Peclet 
number increasing. Contrary, the one-point statistics of 9 
(say, the flatness (6* 4 ) / (9 2 ) 2 ) is independent of the Peclet 
number Pe — L/ra at the limit of large Pe. 

Besides, we find the nontrivial r^-related anoma- 
lous scaling describing the traceless tensor V a 9V p9 — 
d~ 1 8 a p(V6) 2 . We also argue that high-order integrals 
of motion may influence scaling if the pumping is non- 
Gaussian and that the different kinds of L-related anoma- 
lous scaling may be the case starting from the sixth cor- 
relation function. 

The anomalous L-related scaling of the fourth-order 
correlation function of the scalar field confronted with 
the diagrammatic analysis that shows no infrared diver- 
gences at any finite order of a Wyld diagram technique 
p| evidently means that such an anomalous scaling is a 
non-perturbative phenomenon. 



The structure of the article is as follows. We formu- 
late the problem and find the pair correlation function in 
Section U the second section is devoted to L-dependence 
while the third section to r<j dependence of the fourth- 
order correlation functions. The correlation functions 
of the scalar derivatives are also considered in the third 
section. Section || describes possible generalizations and 
Conclusion summarizes the results. 



1. FORMULATION OF THE PROBLEM 

We consider the advection of a passive scalar field 
9(t, r) by an incompressible turbulent flow. The advec- 
tion is governed by the following equation 



(d t - P)6 = 0, 

P(t) = -« a V a + kA, 



(1.1) 
(1.2) 



where both convective and diffusive terms are included 
in the operator P(t), A designates Laplacian here. The 
external velocity u(i, r) and the external source <fi(t, r) 
are random functions of t and r. We regard the statis- 
tics of the velocity and of the source to be independent. 
Therefore all correlation functions of 9 are to be treated 
as averages over both statistics. Averaging over pump- 
ing will be designated by over-line and averaging over 
velocity will be designated by angular brackets. 

A. Basic Relations 

We assume that the source <p is <5-correlated in time 
and spatially correlated on a scale L. The latter means 
e.g. that the pair correlation function of the source 



^i,ri)^(t 2 ,r 2 ) = S(ti - h)x(r X 2), 



(1.3) 



as a function of the argument ri 2 = |ri — r 2 | decays on 
the scale L. The value x(0) = P2 is the production rate 
of 9 2 . Since the pumping is Gaussian then high-order 
correlation functions are determined by x(r). For exam- 
ple, 



0(ii,ri)0(t 2 ,r 2 )^(i3,r 3 )0(i4,r4) 
= x{ri2)x{rzA)5{ti - t 3 )S(t 3 - t 4 ) 
+x{rw)x{rM)S{h - h)5{t 2 - U) 

+xMxW% - h)S(t2 - 1 3 ). 

A formal solution of (|l . l\) is 



(1.4) 



0{t,T) 



dt'U{t,t')(j)(t',r), 



U(t,t') = Tex P(j P(r)dr 

/ dn / dr 2 ... dT n P(n) . . . P{r n ). 

n=0 Jf Jt> Jf 



(1.5) 
(1.6) 
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Since 6 and <fi are related linearly then the pair product 
of the passive scalar field averaged over the statistics of 
pumping is expressed via the pair correlation function of 
the pumping 



0(i,n)0(i,r 2 ) = / dTt/i(t,T)E/' 2 (t,T)x(ria) 

— oo 

t 

= j drT exp(f [Pi W+A(r)]dr) X (n 2 ). (1.7) 



Here and below Ui, Pi designate operators acting on vari- 
ables Yi. Time-ordered exponents in the r.h.s. of (1.7) 



commute with each other (because they hav e di fferent 
space arguments), that allows to rewrite in (1.7) their 
product as a single Texp. 



The next step is to average Texp in (1.7) with respect 
to the statistics of the velocity field u. Following Kraich- 
nan, we consider the case of a velocity delta-correlated 
in time but multi-scale in space. Velocity statistics is 
completely determined by the pair correlation function 



(u Q (ii,r 1 )^(t 2 ,r 2 ))=«5(t 1 -t 2 )y^, 
V aP = V S al3 - /C Q/3 (r 12 ), JC al3 (0) = . 

Here the so-called eddy diffusivity is as follows 



(1.8) 



K *P = Vf r 2 S <x0 _ r a r /3) + D ( d *) 000^-7 . (1.9) 

rf 2 — 7 

where < 7 < 2 is supposed and (. . .) stands for an 
average over the velocity field statistics, isotropy being 
assumed. The representation (1.9) is valid for the scales 
less than velocity infrared cut-off L u , which is supposed 
to be the largest scale of the problem. The marginal "log- 
arithmic" 7 = and "diffusive" 7 = 2 regimes require a 
special care. 



B. Simultaneous pair correlator 

Let us calculate the simultaneous pair correlation func- 
tion of the passive scalar 



/(ri 2 ) = (0(ri)0(r 2 )) . 



(1.10) 



Since we assume the statistics of the velocity field to be 
delta-correlated and Gaussian, then the n-th order corre- 
lator of the velocity field [which appears in the expansion 
of the Texp from the integrand of (1.7)] can be found ex- 
plicitly by the Wick theorem which prescribes to reduce 
any average to a product of pair correlation functions. 
Those terms are summed up into a usual operator expo- 
nent 



Texp [P^t) + P 2 ( T )]dT 



(1.11) 



Vn 



di = V°"(r«)V?V£ + + «j(A, + A,-) , (1.12) 
where = — tj and A$ = Vf. We find from ( fL7| ) 

/(ri 2 ) = J drexp(V-T)£ 12 ^x(ri2). (1.13) 



Integrating the r.h.s. of (1.13) with respect to time (£y 

is a well defined integral operator) we get £12/(^12) = 
— x( r i2) which is an ordinary second-order differential 
equation 



where d r = d/dr 



(1.14a) 



£W = ^y D r i-d dr ^ + r 2- 7 y-i dr f (1 14b) 



previously derived by Kraichnan |15|] . In (1.14t) the dif- 
fusion scale Td was introduced 



2k(2 - 7) 
D{d-1) ' 



(1.15) 



Equation ( 1.14 ) is integrated explicitly. The solution is 
completely determined by two physical boundary condi- 
tions: zero at t — 00 and absence of a singularity at 
r = 0. We consider the pumping correlation function 
close to a step function x( r ) = P2 a± r < L and zero 
otherwise: 



f{r) = P 2 \ g(r')dr', (1.16) 
at r > L g{r) 



• ■ I g(r) = 



2ndl + (r/r d ) 2 -7 ' 

r 1 
2^1 + (r/r d ) 2 ~i ' 



For a general 7, the pair correlation function is expressed 
via a transcendental function 



f(r) = 



P*rl 



d(d-l) \ sin[2vr/(2 - 7)] 



-$[-(r/r d ) 2 - 7 , 1,2/(2-7)] \, 



(1.17) 



- see @, p.29 and Q p.27. In the Richardson- 
Kolmogorov case, d = 3, 7 = 2/3 and one gets 

/(0) - f(r) = y [r 2/3 - ^ /3 arctg [r/r d f^] . (1.18) 



Both (1.17) and (1.18) are valid at any r < L for the 
step-like pumping. At r <ti L, those formulas describe 
the pair correlation function for the case of an arbitrary 
pumping x with the correlation scale L. One can see 
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that the transition from convective to diffusive interval 
is described by a universal function. 

If the Peclet number Pe = L/rd is large, the pumping 
and diffusion scales separate three intervals of scales with 
a different scaling behavior: 



r < r d < L /(r) « P 2 
r d ^r<L f(r) ra P 2 



L-r(2- 7 ) 



7(d-7)(d-l)D 

(1.19a) 

2-7 / L 7 



7(d-l)DVd-7 d 

(1.19b) 

L d (2 - ■y)r 1 ^ d 

^" ^-W-W (L19c) 

Those expressions fully determine the behavior of the 
pair correlation function. 



six at d > 2, which is still too many to hope to find F 
explicitly at all possible distances ranging from zero 
to infinity. Our aim is modest: we are looking for the 
scaling exponents only. First, we shall find the overall 
scaling exponent (4 which describes how T (which is the 
irreducible part of -F1234) scales if all the distances ry 
are multiplied by the same factor. This is the subject 
of this section where we employ 1 j d perturbation theory 
assuming space dimensionality to be large. That will also 
allow us to find the anomalous exponents that describe 
the L-dependence of the correlation function. Second, we 
shall consider the case with one or two distances being 
much less than other ones and find the scaling exponents 
with respect to small and large distances separately, this 
could be done at an arbitrary d. That will allow us to 
fuse some points and to find the r^-dependence of the 
cumulants of second powers of spatial derivatives. This 
is the subject of the next section. 



C. Different-time pair correlator 



The dynamical analog of (1.7) looks as follows (t 2 > ti) 



9(t ll r 1 )6(t 2 ,r 2 ) 

= U 2 (t2,h) f dT#l(tl,T)#2(tl S T)x(ri2) 
J — 00 

- #2(t2,ti)0(*i,ri)0(ii,r 2 ). (1.20) 

Since the velocity field is delta-correlated, then one may 
reduce the average of the product from the r.h.s. of ( 1.20| ) 
to the product of the averages 



{9(t 1 ,T 1 )9(t 2 ,T2))=(U 2 (t2,t 1 ))f(r 13 ) 

= cxp( I t 2 -h\ (V + k)A W12). 



Let us stress that if r<j <C L then Vo ^> k. The point is 
that Vo can be estimated as Vo ~ DLi 2 ^ 1 where L u is the 
scale of the largest vortices. Since L u > L is assumed 
then Vq 3> k. That means that the time dependence of 
correlation functions of the passive scalar 9 is determined 
by the velocity of the largest vortices and is therefore fast. 
Note that in the comoving reference frame this dynamics 
is much more slow flllll . 



2. FOUR-POINT CORRELATION FUNCTION OF 
THE PASSIVE SCALAR 



A. The equation for the simultaneous four-point 
correlation function 



Using the technique developed in the preceding section 
we derive the equation for the simultaneous four-point 
correlation function F1234. Let us form first a four-point 
object averaged over the pumping only by analogy with 
(PI) 



9(t,r 1 )9(t,r 2 )9{t,r 3 )9(t,r 4 ) 

dn I dr 2 I dr 3 / dr 4 Z7i(t, n)U 2 (t, r 2 ) 



f/3(^ / r3)C/4(i,T4)0(Ti,ri)0(r2,r 2 )(/)(T3,r3)0(r4,r4). (2.1) 



Due to (1.4) the object separates onto three parts 

Ql2;34 + Ql3;24 + <9l4;23, where 



irM= / dTdfUiit^Uj&^Ukit^UifafMrijMrki) 



— / drTexp 



dT'iP.+Pj+Pk + Pl) 



Texp df 



T 

xiPi+Pjfj+TexpfJdf'iPk+Pt) 



df 



X(nj)x(rki)- (2.2) 



We begin the analysis of the fourth-order correla- 
tion function ^1234. It will be based upon the equa- 



tion (2.5), previously derived independently by Kraich- 
nan ]20[ ], Sinai and Yakhot pl[ | and Shraiman and Sig- 
gia [ p2| . It is a second-order partial differential equation 
in the space of 4<i variables. Isotropy allows one to di- 
minish the number of variables to five at d = 2 and to 



The commutativity of Pi and Pj at i 7^ j has been taken 
into account at (2.2). 

For the delta-correlated velocity field, one may average 
the time-ordered exponents in the r.h.s. of (2.2) explic- 
itly, first, expanding them into series and, second, per- 
forming the Gaussian decomposition of the velocity cor- 
relators in the same manner as it was done for the pair 
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correlator of the passive scalar in (1.11). That results in 
the following expressions 



(Qi j; kl) = / dre^ c 



df e (- f )^ + e <?-*)£» X (r tf )x(r fc 



(2.3) 



where Cij was defined in (1.12) and C has the following 
form: 



(2.4) 



where Vf = d/drf . The integration over f in ( |2.3| ) gives 

where Aj/(ry) = —x( r ij) was used. Inserting the 
last expressions into (2.3) and collecting all the terms 
(Qij;ki) we find for the full four-point object -F1234 = 
(0(ri)0(r 2 )0(r 3 )0(r 4 )): 

^1234= / dr exp[(t-r)£] f(nj)x(rkl), 

where the notation {r/} stands for a pair of non- 
coinciding (i 7^ j) site indices. That leads to the following 
equation 

-£^1234 = /(ri 2 )x(r 34 ) + /(r 34 )x(r 12 ) + /(n 3 )x(r 24 ) 
+/(rs4)x(ri3) + /(r-i4)x(r- 23 ) + /(r 23 )x(n 4 ), (2.5) 

with the r.h.s. expre ssed in ter ms of the functions f(rij) 
already found: see ( 1.16 1.19 ). Note that at a conflu- 



ence of the points — > r 7 the r.h.s. of ( |2 . 5|) tends 



to a constant. Using (1.14) one can obtain an equa- 
tion for the irreducible four-point correlator 1^1234 = 
((0(ri)0(r 2 )0(r 3 )0(r 4 )» = i^ 234 -/i 2 / 34 -/ 13 / 24 -/i 4 / 23 : 



234 



12:34 



13;24 



14;23, 



(2.6a) 



ij;kl 



JC aP (v a ) - JC ap {y lk ) + JC aP {v ]k ) - JC a0 (r 3l ) 



xVg.Vg,/(ry)/(r w ), 



(2.6b) 



where Vf. 



The operator £ ( ^.4| ) is negatively defined. To prove 
this we represent £ — nCdi / + C c where the convective 
part could be written as £ c = ( [VJ u(rj)Vi] 2 ) u - Here we 
designate (. . .)„ as an average over the velocity field u(r) 
random in space - from the viewpoint of initial u(r, t), 
the averaging is over the instant configuration. Incom- 
pressibility guarantees that £ c is Hermitian and nega- 
tively defined as well as Cdif- Therefore, £ has a con- 
tinuous non-positive spectrum and the density of states 



regular at zero. The last statement follows from the in- 
equality \\n£dif + C c \\ > || KjCdi/ 1| and the absence of sin- 
gularity in the density of states for Cdif- We can thus 
conclude that the equation — CT = $ is well defined for 
r and $ that do not grow at infinity. 

The operator C is scale invariant if all the distances 
Tij are either much larger than or much smaller. The 
right-hand side is scale invariant if all the distances are 
either larger or smaller than L. We thus could divide 
our space of into three domains where the scale in- 
variance of r(ry) is to be expected. It is natural to ask 
now a simple question: what prevents us from making 
the statement that the scaling exponent of the solution 
is equal to the difference between the exponents of the 
right-hand side and of the operator CI Of course, the 
solutions with such "naive" scalings exist in all three re- 
gions (we call them forced solutions Ft orc ). The problem 
is that to match those particular solutions at pumping 
and diffusion scale, it may be necessary to include into 
the full solution the zero modes of the operator: 



F — Ff orc + Z 



(2.7) 



Here the zero mode Z may have a scaling different from 
that of Ff orc . To avoid misunderstanding, note that 
since the operator C is non-positive then it cannot have 
a global zero mode that satisfies boundary conditions. 
The parts of the solution with an anomalous scaling may 
nevertheless be considered as zero modes within separate 
domains. At r.y <C L, for instance, one should worry 
about zero modes that appear due to matching condi- 
tions at ~ L, such modes are allowed to grow with 

. At an oversimplified level, this is what happens when 
we ob tain the con stant L-dependent term in the expres- 
sions ( 1.19a , 1.19b ) for the pair correlation function. In a 
multidimensional space, operators may have zero modes 
much more complicated than a constant and, indeed, the 
operator C does have an infinite number of zero modes. 

We cannot yet find the zero modes and solve the 
matching problem analytically at arbitrary d. Fortu- 
nately, for d > 2 (generally, for d > n — 2 where n is the 
number of points in the correlation function) we found 
the representation of C that allows one to represent Tfcj) 
(with all distances of the same order or the main term 
if some distances are smaller than others) as a power se- 
ries with respect to some numerical parameter (1/d in 
the convective interval). The coefficients in the series are 
functions of all having the same scaling exponents and 
logarithms that appear from expanding anomalous expo- 
nents in powers of 1/d. By analyzing those functions at 
the next subsection, we establish the overall scaling prop- 
erties of T(rij) in the convective and diffusion intervals. 
The possibility to sum the logarithms into a power func- 
tions is provided by an explicit scale invariance of the 
equation that determines the zero mode in the convec- 
tive interval. The consideration of the subleading terms 
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at the case of some distances being small in comparison 
with other ones requires a special approach that will be 
developed in Section ||. 



B. Representation of Tetrahedron's Lengths 

To establish overall scaling properties of the irreducible 
fourth-order correlator T it will be convenient for us to 
use a special representation of the equations fl2.5|]2.6| ). In 
an isotropic case, the four-point correlation function T is 
a function of six distances rtj between the points, d > 2 
being assumed. The case d — 2 where there are only five 
independent variables deserves special consideration. In 
the variables r, 



| r i,- 1 , the operator £ is a sum of two 



parts, C = Cq + C\\ 



C = 



ilVr 1 ^ (r 2 ~ 7 +r 2 ~~ f )r d - 1 



i>j 



M 2(2 - 7 ) ^ m 13 3n> r jn dr t3 dr 3n 



1-7 



d 2 



l 



D 



4 FijTimTjn 



^•^^-^^2/2,2 2 2 \ 

2 ,y 13 V in 1 jm ij ran) 



,1/2 1 2,1-1 'J _> , 

a — (r ■ + r ■ — r )[r • ■ + r • — r- 

1 %j 1 im jva/Vij 1 jn in/ 

r>-»2 J !Y*^ /v*- 

'11 I ' 1 til ' • 



-K 



2 — 1 — j.2 j,2 



jm/v ij 1 j7 

d 2 



d 2 



(2.8) 



Here the summation is performed over subscripts satis- 
fying the conditions i ^ j and m ^ i, j , n 5^ 

What is nice about that representation is that if we 
omit C\ from £ then the reducible part Ff^ = /12/34 + 
j 13/24 + /14/23 o f th e fourth-order correlator appears to 
be a solution of (|2.5|). It is 



^ (0) = f( r v)f( r ki) = (in the conv. int.) 
3(2- 7 ) 2 L 2 ^ , (2- 7 ) 2 ~ 



7 2 (d-l) 2 (d- 7 ) 2 ' 7 2 (d-l) 2 d 2 +F/on 



2(d + 7 ) v lJ 



fcZ 7' 



(2.9) 



F 



(2-7) 2 E 



.27 



fore - 



(2 - 7 ) 2 £T V rX 

- " (2.10) 



2 7 2 (d + 7 )(d- l) 2 d 7 2 (d- l) 2 d(d- 7 ) 



which is the zero approximation of (2.7). Here -Zj^/y 
is the zero mode of £0 while that is remarkable that 
a direct check shows Ff orc to be a partial solution in 
the convective interval of the full equation —CFf orc = 
/( r i2)x(^34) + • • ■• It means, particularly, that if we de- 
velop the iteration procedure in £q 1 C\ then all the terms 
appeared at higher steps will enter zero modes of the full 
operator C but not the forced term found. Another re- 
markable feature of the zero step is an absence among 



the terms three-point zero modes of Cq, like Zij^j.. They 
will appear on the next step only. 

The basic fact is that Cq oc d 2 while C\ cx d as d — > 00. 
Assuming 1/d to be a formal small parameter, we shall 
implement the iteration procedure with respect to C\ 



£ r (1) = -c x f^ , £ r (n) = -£ir( r 



(2.11) 



This procedure leads to a representation of T as a series 
over the powers of 1 /d: 



r = J2(-Co 1 c 1 ) n F^ 



(2.12) 



which is actually the series for the nongaussian part of 
the zero mode. 

One may note that the four-point correlation function 
is defined not in the whole six-dimensional space of 
but rather in the physical subspace restricted by triangle 



inequalities 



> r^. That makes no additional 



difficulties since the solution we shall find satisfies all 
boundary conditions in the physical subspace. In addi- 
tion, the answer is expressed in terms of powers of £ ( j" 1 £i 
which do not have singularities at the boundary of the 
subspace. 

At each step of the iteration procedure we should solve 
the equation of the (2.11) type. Such equation is much 
easier to analyze than e.g. (2.6a). The point is that the 



operator £0 is a sum of the six sub-operators C^ from 
( 1.14b| ) each depending on the single variable only. 
This particular form of £0 enables us to analyze the zero 
modes and establish the necessary properties of its resol- 
vent. The solution of the equation (2.11) is expressed via 
the corresponding integral kernel as follows 

rW(f) = j dtJlj ' dr^niW^^ir 1 ), (2-13) 

i>j 



where r designates the set of six variables and the 
integration is performed over time t and six separations 
r\; . The resolvent 1Z can be represented as a product 



(2.14) 



i>3 



where the function TZ^ P ' (t;r,r') satisfies the equation 

(d t - Ci p) )TZ^(t;r,r') = S(t)5(r - r') , (2.15) 

with the condition TZ^(t;r,r') = at t < 0. Note the 
essential property of 1Z simplifying the subsequent anal- 
ysis: it is independent of the pumping scale L. That 
means that there is only one characteristic length in 1Z: 
the diffusion scale introduced by (1.15). 

Now we are going to establish the properties of the re- 
solvent 1ZS p \t]r,r'). Let us remind that there are two 
terms in the operator C^ (1.14b) which are of the same 



G 



order at the diffusion scale rd- Solving (2.15) at the dif- 
fusion and convection intervals, we get for the resolvent 



at r, r', Vtn <C r c i 



(2.16a) 



n (p) 



1 



exp 



at r y,(Dt) 1/ri > r d 



d/2-l 



rr 



r X-d/2 r ld/2 



exp 



(2-7) 



j(d-l)W 
(ri + r n ){2- 1 ) 
-y 2 (d-l)tD 



r (7-<0/2/(d+7-2)/2 x 



d/j-i 



(2.16b) 

/ 2(2- 7 )(rr') 7 / 2 ' 
j 2 (d-l)tD 



We shall need also the asymptotics in the mixed limit 
r 7 r d r', (Dt) 1 ^ , for r and r' lying in the different in- 
tervals (diffusion and convective ones respectively) . This 
asymptotics has to match with one following from (2.16t ) 
in the sub- interval <C r <C r', (Dt ) 1 ' 1 . It gives an idea: 
to improve the expansion of ( [2.16b| ) with respect to r/r' 
introducing a function y(r): 



at 



r, r d «r',(Z?i) 1 ^ 



n {p) 



(2-7) 
j 2 (d~l)tD 



ld/l Jr >(d-X) 



exp 



r'7( 2 - 7 ) 



-1 



2(2 - 7)r' 7 
- 1)7* 



r(d/ 7 ) 

2(2- 7 )y(r) 



7 2 (d- 



1)£D 
(2.17a) 



In the limit <C r, the function y(r) should pass to 
r 7 as it follows from (2.16b). Substituting (2.17a) into 



(2.15) and solving the resulting equation for y one gets 
an expression with one unknown parameter: 



y(r)=y{0) + (d-l)D 



xdx 



D(d- l)x 2 -~« + 2n{2 ~ ~f) ' 
(2.17b) 



Generally, the integration of the second order differential 
equation produces two parameters, but here one of them 
has been already fixed by the condition of the finiteness 
of y(r) at r — > 0. The only scale which can determine the 
dimensional parameter y(0) is r^. It gives the following 
estimate y(0) ~ rj and, thus, closes our analysis of the 
mixed asymptotics of the resolvent at r, <C r', (Dt) 1 ^ 1 . 



C. Overall Scaling of the Fourth-Order Correlator 



The analysis of the series (2.12) enables us to estab- 
lish the scaling behavior of the irreducible part of the 
fourth-order correlator T in different regions of scales. 
The crucial point is that all terms of the series ( 2.12] ) have 
the same scaling behavior as I^ 1 ) up to some logarithmic 
functions that are our main concern. Indeed, we should 
worry if some logarithms could appear at each step of 



the iteration procedure (2.11), summation of powers of 
logarithms can produce anomalous exponents changing 
the index of the whole sum in comparison with the first 



term of J2.12Q . 

We shall analyze in this Section only logarithms 
ln(L/rij). Besides, the logarithms of the ratios of separa- 
tions ln(rjj/Vfe;) may arise at any step of the iteration pro- 



cedure (2.11) which could change the behavior at small 
ratios. This happens only in subleading terms where the 
coefficient at logarithms are proportional to a power of 
the small ratio. One may worry, nevertheless, if the log- 
arithms could be summed into the large exponent which 
compensates the small factors. That means that the case 
of small ratios needs a separate non-perturbative analy- 
sis which is done in Section ||. The results of the analysis 
show that after re-summation the subleading terms re- 
main the subleading ones. Therefore, the logarithms of 
the ratios can be neglected at the general investigation. 
Below in this section, we will speak about the overall 
scaling of the zero modes that is about their exponents 
in terms of the ratio L/r. 

First, we show that the zero modes of Cq cannot con- 
tain logarithms ln(L/ry ■). It can be proved using an "an- 
gular" representation in -space introduced as 



r 7/2 
'12 



R l/2 



Hi 



r 7/2 
13 



£7/2 



n 2 . 



(2.18) 



where R is a "radial" variable and n is a unit 6- 
component vector. In terms of those variables the op- 
erator Cq at R ^> rd is written as 



(d-l)D 
(2 - 7)i? 7 



R 2 d 2 R + (6d + 1 - "f)Rd R + T 



where T is an angular operator: 

- 7(1 -7 + d) / d 

T = ~ 2^ ( n »> ~ 6rim >- 



(7/2) 2 £ . 

771—1 

d 



2/7- 



m— 1 

dn m 



dn r , 



+ 



,1-2/7. 



dn r - 



(2.19) 



dn r . 



dn r . 



6 a 
> n m m — 



The consequence of ( 2.19 ) is a possibility to look for the 
zero modes in the following form: 



Z n = R a "Y n {n), 
TY n (n) = X n Y n (n), 
a 2 + (fid - 7)a n + A r , 



0. 



(2.20) 
(2.21) 
(2.22) 



The roots of the equation (2.22) determine exponents 
in the expression (2.20) for the zero modes which ap- 
pear to be power functions [like fl2.9|)]. Dangerous loga- 
rithms would occur if the roots of (2.22) coincide. That 
corresponds to the values a rhCr — — (fid — 7)/2, A„ 



7 



(6tf — 7) 2 /4, which is impossible since all eigenvalues A„ 
of the operator T are negative. The property follows 
from the fact that all eigenvalues of Cq are negative and 
Cq is reduced to T for i?-independent functions. Non- 
degeneracy of the eigenvalues guarantees the absence of 
the logarithms ln(i/ry) in the zero modes of Cq. Note 
that the coincidence of a n for two different A„ (i.e. the 
existence of two zero modes with the same scaling yet 
different angular structures) apparently does not lead to 
the appearance of logarithms as well. 

To establish the exponent of the overall scaling, let us 
analyze the integral (2.13). First, we show that at any 
step of the iteration the integral converges at r' ^> L so 
that this region gives negligible contribution if we con- 
sider separations r.y < L. To do that, one should find 
the behavior of 1^") at large scales. The decay of $ 



from (2.6b) is different for different direction in the six- 
dimensional space {rij}. By analyzing the most slowly 
decay ing case we get $ ~ B r 1 w here R is determined by 
( 2.18 ). Now we can analyze ( 2.13| ) and using the resolvent 



which is determined for all scales one can directly show 
that the region R' ^> L gives a negligible contribution. 
Then the analysis of in the region R^> L shows that 
rW decays like $. That means that the above scheme is 
reproduced at any step of the iteration procedure. 

Consider the case where all rv, are in the convective 



interval. Then we can use (pU6b) for The contri- 



bution to associated with r' < r is L-independent 



One can check using (2.161) and (1.19b) that the integral 
in (2.13) converges at small r' already in the convective 
interval. That means that the contribution associated 
with r' from the convective interval is also r^-independent 
and it is thus given by the simple dimensional estimate: 
P 2 / D 2 multiplied by a scale invariant function of with 
the exponent 27. That contribution could be neglected in 
comparison with the terms proportional to positive pow- 
ers of L which will b e found below. Using the expressions 
analogous to fl2.17| ) for small r' one can also check that 
the contribution to the integral due to the diffusion re- 



gion is finite and negligible. The main contribution to 
r^") at the integration over r'^ in (2.13) is given by the 
region ~ L, the characteristic time being estimated 
as t ~ L^/D . It follows then from fl2.16b| ) that for 
r <C L the corresponding expression for "R~ffi is a regu- 
lar expansion in r 1 , what means that the contribution 
to associated with the considered region r\- ~ L is 
also a regular expansion in rJ-jL 1 . Only the first terms, 
proportional to positive powers of L, are important, the 
subsequent terms proportional to negative powers of L 
can be neglected in the convective interval. 

The zero term in r 7 -expansion of ( |2.13 ) is a constant 
which can be estimated as P 2 L 21 /D 2 . The first term 
(that would be proportional to L 7 r 7 ) is actually equal to 
zero. That could be understood as follows. As it is seen 
from fll.191: , 2.8) in the convective interval, the r.h.s of 
( |2.11 ) docs not depend on L: (1.19b) contains only an 
L-dependent constant, these constants has to be killed in 
( |2.13 ) since the operator (2.8) contains only cross deriva- 
tives. That L-independence of the r.h.s. means that an 
L-dependent term in rW should be treated as a zero 
mode (up to a possible logarithmic factor) of the oper- 
ator Cq (in the convective interval only!). Zero modes 
with the scaling r 7 exist (e.g. rj 2 — r 3A ) but one can- 
not construct such mode symmetric in permutations of 
the points 1,2,3,4. Thus the presence of such zero modes 
would violate the intrinsic 
lem. 

The next term behaves o 



<-> Tj symmetry of the prob- 



r27 



If contains loga- 
rithmic factors multiplied by r 27 then the normal scaling 
is renormalized into an anomalous one. Let us consider 
the general expression for r' 1 ) 



dt / d^R^r,?)^^^) , (2.23) 



r«(r) 



expand it in the series in r 7 and calculate (for instance) 



the terms proportional to ^ 



7 7 

12' 34' r 12 r 13 



and 



, 27 



tively: 



Ar (1) - / 
12,34 - ' 



IA\ jAx(tD/V), AT^ 13 = rJ 2 rJ 3 J jA^lD/V), Ar« 12 = r\l J ^A 3 (tD/L^), 



A 1 {t/L"<) = dr' 12 n^(t;0,r' 12 ) dr> 34 K^(t;0,r> 3i ) dr' 13 K^(t; 0, r[ 3 ) / dr' 14 7^> (i; 0, r' 14 ) x 



x dr' 23 n^(t;0,r 23 ) dr' 2 ^\t-,Qy 2A )C^)F^{f), 



A 2 (t/r)= dr' 12 n^(t;0,r' 12 ) dr' 13 TZ^(t;0,r' 13 ) dr' 3i U W (t; 0, r' M ) / dr' u 1l W (t; 0, r 14 ) X 



x dr' 23 n^(t;0,r 23 ) dr 24 TZ^(t; 0, r^C^F^i?), 



~(p) 



A 3 (t/r)= dr' 12 K (t;0,r' 12 ) dr' 34 K^ (t;0,r' 34 ) dr' 13 K « (i; 0, r' 13 ) / dr' u TZ^ (t; 0, / u ) x 



x dr' 23 1lW(t-Ar 23 ) dr' 24 Tl^(t;Oy u )M^F i0) (^ 



respec- 



(2.24) 



(2.25) 



(2.26) 



(2.27) 
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K {p) {t;0,x) 



2-1 
7 2 (d- l)iL> 



d/j 



T[d/ 7 ] 



exp 



(d- 1)Y 

(2-7) 2 
2(d- 1) 2 7 4 



_ 7 2-7 

1 7 2 (d-l)iD 
(2-7)xT 



~(p) 



^"'(t;0,a:) = ^ (p) (t;0,a;) 



d(d - l)jtD 
2(2 - 7 )a; 7 



(2 - 7 ) 2 :r 27 



d(d-l) 7 iL> d{d- l) 2 (d + 7) 7 2 t 2 D 2 _ 



(2.28) 
(2.29) 
(2.30) 



While integrating over time at tD <C one can 
expand the functions Ai{tD / L 1 ) in a regular series in 
tD/LP: 



Ai(tD/L 7 ) = C { 1) + C^tD/IP + 



(2.31) 



The crucial point is the presence of the zero term in the 
expansion. Were Cq be equal to zero we would imme- 
diately figure out the part of the well known zero-mode 
from the divergent integral over t (~ Cf 1 ) in ( p. 24 ), the 
respective term in (2.24) being L- independent. As we 



show below, C$ ^ so that one gets the logarithm 
ln[L/r] in the respective (~ Cq) contribution to AIj . 
This term stems from the forced solution of the equation 
CoT^ = — CiF^ (it is proved already that there are 
no logarithms in zero modes of Cq). From the viewpoint 
of the initial equation (2.5), those logarithms stem from 
the zero modes of the full operator C — see ( [2. 10 ) and 
below. 

(i) 

Let us calculate Cq directly. First of all one finds 
from (2 



■AF (0) =- 

(2-7) 2 
" 8(d- l) 2 d 2 



2-7 



7 (d 

E 



l)d 
.d4 



7 „7 
2 '34 



7 7 
r 13 r 24 



7 2 / 2 

— r ■ r 



2 2 
T • — T- 



,/2i2 2 \ / 2 , 2 2 \ 



—7 7 — 
ij im 



2 7-' 

r ■ 

jn 



(2.32) 



where the summation is performed over all the sets of 
different subscripts {i,j,m,n} = {1,2,3,4}. F^°\r) and 
respectively the integrands in ( 2.28||2.3C| ) are sums of a 
huge number of very simple terms that have been cal- 
culated by Mathematica. We do not present here the 
cumbersome expressions obtained formally at arbitrary 
7 and d since the physical meaning could be ascribed to 
them only when £i/£o-expansion is a small-parameter 
expansion. The below formulas show that there are two 
cases where this is so: 
d — > oo: 



Cq go to zero at 7 — > 2 and at 



at d — ► 00 Cn 



at d — > 00 C, 



(2) 



at d — > 00 C, 



(3) 



4(2 - 7 ) 3 
7 3 d 5 

2(2 - 7 ) 3 
7 3 d 5 

2(2- 7 ) 3 ; 
7 3 d 5 



(2.33) 
(2.34) 
(2.35) 



at 7^2" 
at 7 -> 2" 
at 7 -> 2" 



C, 



c, 



(1) _ (5d-2-d 2 )(2- 7 ) 3 
2(d-2)(d- l) 3 d 3 ' 

(2) _ (2 + d 2 -4d)(2~ 7 ) 3 
^ 4(d- 2)(d- l) 3 d 3 ' 

(3) ^ (2-7) 3 
4d 2 (d- l)(d-2)' 



(2.36) 
(2.37) 
(2.38) 



It is natural that our £q~ £1 -expansion gives small cor- 
rection in the limit d — ► 00 while we should admit 
that the same behavior in the limit 7 — > 2 is surpris- 
ing for us. That means that /Cq -1 /^ -expansion might 
be meaningful also when d is arbitrary while 2 — 7 is 
small. To establish that, one should check whether non- 
logarithmic terms in are proportional to (2 — 7) 3 
as well, which is beyond the scope of the present ap- 
proach. A nontrivial task to develop a consistent (2 — 7)- 
expansion will be the subject of further publications. As 
7^2 velocity spectrum diverges at small scales, one 



should introduce ultraviolet cut-off r v in ( |l.9[) so that 



r 7 /(2 — 7) — * \n(r/r„). The scalar correlation func- 
tions contains logarithms already at 7 = 2 (in the paper 
|i~0f , the limit 7 — > 2,D/(2 — 7) ^const has been consid- 
ered where such logarithms are absent). 

Contrary the perturbation theory in 1 / d which is the 
main subject of this Section is regular and uniform with 
respect to distances r^ . Note that we assume 1/d to 
be the smallest parameter in the problem, in particular, 
1/d <C 7 so that our results do not describe the limit of 
small 7 which will be the subject of further publications. 

We have found all the logarithmic terms entering the 
fourth-order correlator on the first step of the iteration 
procedure: 



f 12 ,34=^ 2 °; 34 +7 K ] £ 



' ij ' kl 



ln[L/r 



"Cf £ r l r l ln ^/V] + C< 3) ]T r% HL/r] 
-L-independent or non-logarithmic terms + 



(2.39) 



where the summations are performed over all the possi- 
ble combinations supposing all the indexes k, I to be 
different. The expression in front of the logarithms in 
( p. 39 ) is a sum of the well-known zero modes (~ r 2j ), 
that is guaranteed by the identity that could be directly 
checked 
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<#) + _£_(cf+4Cf)=0. (2.40) 

It agrees with the absence of any logarithms among zero 
modes of £ . 

We have thus found the first (in 1/d) term of the ex- 
pansion of the zero-mode part of the solution of the full 
operator 



Z = 



(2-7) 5 



7 2 d 4 



E z <.^ ■ - (<w e^- i ^ L /^ 



(2.41) 



Note that there are two types (Zij^ki and Zij^k) of zero 
modes of Co possessing scaling 27, only one of them 
(Zij t ki) is present on the zero step of the iteration proce- 
dure. It is clear that ( 2.41 ) cannot be described in terms 
of a single anomalous exponents A(d) because the expres- 
sion in the last brackets does not coincide with ^ Zij t ki- 
We thus come to the conclusion that Z should consist 
of at least two terms with different scaling. On another 
language it means that one gets a degenerate case. In- 
deed, The full operator £, the bare operator C and the 
perturbative one C\ are scale invariant. The bare zero 
modes scales similarly. And it is very essential that the 
first perturbation step should dismiss the degeneracy: on 
the first step of the iteration procedure it should appear 
two zero modes 



Z2 = E 



(aZ ijt ki + fiZijfk) (l - Ai ln[r/i]) 



(2.42) 



((1 - a)Z l]M - f3Z ijtik ) (1 - A 2 ln[r/L]) 



(2.43) 



the logarithmic terms appeared are proportional to the 
bare zero mode ^ Zi n j n only. It is a manifestation of 
the fact that the functions on a triangle constitute an 
invariant subspace of the full operator C: L acting on 
an arbitrary function of three distances from a trian- 
gle ri n ,rij,rj n produces a function of the same three 
distances again. Thus the cancelation of a coefficient 
before the logarithm proportional to another four-point 
zero mode will occurs on all the higher steps of the it- 
eration of ^ Zi n j n too. Thus, the resulting zero mode 
(say Z\ with a = 0) is scale invariant with the following 
asymptotics of the anomalous exponent at d — * 00 



Ai 



(2- 7 )(2 + 7 ) 
2d 



(2.46) 



The asymptotics of all the rest coefficients defining the 
second (mixed) scale-invariant zero mode Z2 is restored 
from (|2?4ll,^42|,^"43|j2.44^46|) 



10 + 7' 



(2.47) 



Note that the above £j~ 1 £o iteration procedure was not, 
strictly speaking, a direct l/<i-expansion: not only ma- 
jor terms ~ d 2 but also subleading ones ~ d have been 
included into Cq. That made it possible to have both 
unperturbed operator and perturbation of the same (sec- 
ond) order. To reinforce the results and obtain the ex- 
plicit solution for Ti (not only it's logarithmic part as 
above), let us briefly describe the direct 1/d-procedure 
where the main part of C in the convective interval is the 
differential operator of the first order 



C' 



D 



(2.48) 



i>j 



such that 



(2~7) S 
7 2 d 4 



(Z1 + Z2), 



(2.44) 



where Z was found above in ( 2.41 ). 

To extract an additional information that will fix all 
the numbers a, (3, Aj we should calculate first-order log- 
arithmic correction — C^ 1 C\ ^2 %in,jn to the second bare 
zero mode ^ Zm,jn ■ First of all one finds 



(rf-l)7 2 
2(2-7) 



E 1 



2 2 
r- — r 

ij in 



2 \ 7-2 
in)' in 



E 



7-2 7-2 

in jn 



^ + 1 72/2 1 2 21 



1 I 2 1 2 2 \t 1 , 2 2 \ 

q V i>3 m jnj \ %j 1 jn in J 



(2.45) 



yet we will directly account for the necessary boundary 
conditions while integrating over characteristics below. 
The zero modes of ( p.4£ ) which we are going to iterate are 
d — > 00 limits of J2 Z i]M and Z l]iik : Z — J] [2r7-r^ - 



V 1.1 ' 



the first-order corrections to the bare zero modes one has 



E[2 



7 7 



„ 2 7 



*)]• 



To find 



to solve the differential equation 



A = ^{E(^ 



1 



' ' in/ 



1-7 



7 dry) 



d 2 



2 £^Vm ij jnJ Br dr 

1 r 2 ~ 7 fl 2 

/ 2 2 _ r 2._ r 2 \ O 

A Z^i r . r . \'m~'jm 13 'ran) a a 
^ ' m ' in L// im u 1 1 



(2.49) 



(2.50) 



Performing calculations analogous to what was done at 
the calculations of the coefficients Cg one finds that all 



where £[ is the part of the full operator (in the convec- 
tive interval) proportional to 0?, note that it stems both 
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from Co and C\. The equations ( 2.49| ) are integrated by where .By should be considered as constants, and Z[ it- 
characteristics, for example self is defined upto zero mode s of the bare operator C' . 

The result of integration (2.51) at an arbitrary 7 is bulky, 



Z' 



2=7 



t + Bi 



(2.51) 



we present here the explicit expression in the simplest 
possible case of 7 = 1 



l>3 



ln[L/r 



im ' jn 



{Tjn Tij ) [ (^jn Tin) ^Tjn ^jm) ~h (^jn ~l~ ^*mn) ] 



ln[L/r, 



^(7 = 1) 



" d ^ r ' lJ d ^ 



i>3 



rik - r kj 



(r%j + r k j - 2r ik ) \n[L/r lk ] - (ry + r ik - 2r kj ) ]n[L/r k j] 



One can directly check that those are the solutions of 
(2.6) at the first order in 1/d. Being interested only in 



finding the overall scaling we can just put all the sepa- 
rations under the logarithm to be the same. Generally, 
one can extract the logarithmic parts of Z[ , Z[ directly 
from the r.h.s. of ( p. 51 ) putting formally the upper limit 
of the integral to be very large and looking for the loga- 
rithmically divergent terms 

Z[ = ^p^(2Z' - Z' Q ) ln^^+non-log. terms, (2.52) 

-t 2 - 4 - { L\ 
Z[ = — Z Q ml - j +non-log. terms. (2.53) 



This is t he direct analog of ( 2.41 ). The formulas 
( [2.52 2.52 ) imme diately give the leading d — > 00 asymp- 
totics ( 2.46 , 2.47 ) obtained above from the general itera- 
tion procedure. 



The str ucture of the resolvent ( 2.16b ) at small r deter- 
mined by ( 2.17 ) shows that for there exists the limit 
k — ► 0. This property being reproduced at any step of 
the iteration procedure enables one to construct a solu- 
tion r( co ) which is the limit of T at k — > 0. The function 
p(co) jg c \ ose t y in the convective interval and remains 
finite at any ry — > 0. The behavior of r( co ) is deter- 
mined by the estimation (2.54) where r is the maximal 
value among nj . 

We thus come to the conclusion that for separations r 
from the convective interval r<j <C T <C L 

r(r) - T(0) ~ P 2 V 7 ( J L/r) A /.D 2 , T(0) ~ P 2 L 27 /L> 2 

(2.54) 

up to dimensionless constants depending on d and 7 and 
the anomalous dimension A is given by the largest expo- 
nent (A2) which is asymptotically 



A = A 2 = 4 



2- 7 



0(l/d 2 ). 



(2.55) 



Finally, we treat the case where all separations are in 
the diffusion interval. The contribution to associated 



with r' ^> rd gives the constant ~ L 21 '(L / V d) A P% I D 2 ■ 

rd has to 



The contribution associated with scales 



< 



be analyzed carefully since for all separations from the 
diffusion interval the typical integration time t in (2.13) 
is also characteristic of the diffusion region. We use 
fl2.16a| ) for and ( |1.19cj ) for /. The estimate for the 
r-independent contribution gives ~ r 2 ? P 2 /D 2 . To ana- 
lyze a r-dependent contribution associated with r' 



< 



it is convenient to return to the initial form (2.6) of the 
equation for T. At the analysis we will not divide C into 
Co and Ci, the consequent conclusions are consequently 
non-perturbative. In the diffusion region, C is reduced 
to the sum of Laplacians. The r.h.s. of ( [2.6b| ) is pro- 
portional to r 4-7 , that generates the forced part of the 
solution proportional to r 6 ~ 7 . To find the principal r- 
dependent term one has to compare the forced solution 
with zero modes of the sum of Laplacians. All the modes 
are well known: they are constructed from even powers 
of Tij . There are no modes with second powers of r y in T 
since they are not invariant under permutations of points 
17 The zero mode of C of the fourth power possessing 
the permutation symmetry: 



E 



2 2 



lJ nm 2{d+2) y « 



(r 4 + r 4 ) 

vij 1 nmJ 



(2.56) 



turns out to produce larger contribution than the forced 
term. Thus, we conclude that the principal term in T 
is of the forth power in r. The parametric dependence 
of the coefficient at the power can be easily established 
from a matching at r ~ Td' 

at r <C r d 

T - (dL 27 + C'^Pi/D 2 ~ r 4 r 2 r\LlrdtPllD 2 . 

(2.57) 

The consideration of this Section shows that the de- 
gree of non-Gaussianity increases as the distances de- 
crease in the convective interval and eventually comes 
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to a constant when the distances are in the diffusion in- 
terval. For example, the flatness factor is proportional to 
(L/r) A in the convective interval L ^> r ^> and it is 
an r-independent constant of the order of (L/rd) A in the 
diffusive interval r c i 3> r. 



3. FOUR-POINT OBJECTS WITH STRONGLY 
DIFFERENT SEPARATIONS 

In Section || we established the overall scaling behavior 
of the four-point irreducible correlation function 1^1234. 
In particular, one may check that at any step of the it- 
eration procedure ( 2.11 ) the contributions to 1^1234 as- 
sociated with possible small values of the separation ra- 
tios are negligible in comparison with the leading term. 
To show that the summation over n could not convert 
subleading terms into leading ones, here we check the 
self-consistency of our approach nonpertubatively in 1 /d 
using the full operator C. In the next subsection, we show 
that if any separation p is much less than other ones then 
the p-dependent contribution to Ti234 is much less than 
the leading term. This property enables one to find such 
two-point objects as ((6*i — #2) 4 ) which can be extracted 
from Ti234 by fusing some points. Besides, to establish 
the properties of the correlation functions ((ei#§)) and 
((£162)) one should know not only the limit of Ti234 for 
coinciding points but also the dependence of 1^1234 on 
distances between nearby points in groups strongly sep- 
arated from each other. It is the aim of this section to 
obtain the dependence. Here, we shall use perturbation 
theory with respect to the ratio between small and large 
distances. Calculating the corresponding contributions 
to Ti234 one can obtain the two-point correlation func- 
tions of products of different spatial derivatives. We shall 
demonstrate how an absence of an r^-related anomalous 
scaling in 1^1234 prescribes it for ((ei^)) an d (( e i £ 2))- 



A. Four-Point Correlator with a small Separation 

Suppose that among separations ry there is the sepa- 
ration p, say r34, which is much smaller than other sep- 



arations. Then the main term in the operator C (2.4) is 
associated with derivatives over the vector p = T3 — T4. 
This term can be written as 



(3.1) 



where V" = d/dp a . To solve the equation ( |2.6a| ) for 
Ti234 in this case, we can formulate an iteration proce- 
dure where C p is treated as the main contribution to C. 



Namely a solution of ( 2.6a ) can be represented as 



1234 



r (ri 3 , r 23 ) + ^2 r «( r i3> r 23, p) , 



(3.2a) 



-£ p r n = $ n (r 13 ,r 23 ,p) (3.2b) 

$1 = CT (rx, r 2 ) + $i 2;3 4 + $ 13:24 + $14;23i (3.2c) 

$„ = (£-£ p ) r («-i), n>l, (3.2d) 

where ri3 and r23 are separations much larger than p 
and To is p-i nde pendent part of T. The procedure in- 
troduced by (3^2) can be considered as a series over the 
small parameters P/V13, p/r23- We impose two boundary 
conditions on r„: to be of order of ro at p ~ 7-13, r 2 3 and 
to remain finite at p — 0. 

The r.h.s. of (3.2b ) can be assumed to be known. A 
solution of (3.2b) satisfying the imposed boundary con- 
ditions has the following form 



r„(ri,r 2 ,p)= J dt J dr7£(t;p,r)$„(ri,r 2 ,r) . 



(3.3) 



where 1Z is the kernel of the resolvent of the operator C p : 

(d t - £ P )K(t; p, r) - 5(t)S(r - p) . (3.4) 

Since the operator C p explicitly depends on p the resol- 
vent 7Z is not a function of the difference | r — p | only 
as it would be in a homogeneous case. The symmetry 
admits to treat TZ as a function of p, r and x = pr/pr. 
It is convenient to expand the function in the series over 
Jacobi polynomials [p3[ which are the eigenfunctions of 
the ^-dependent part of C p : 



n{p 1 r 1 x)=Y j ^ 2n \p 1 r)P^\x) 



(3.5) 



where v = (d — 3)/2 and only the polynomials with even 
numbers participate since TZ is an even function of x (at 
d = 3, Jacobi polynomials turn into Legendre ones). One 
can obtain from (3.4) separate equations for TZ ty2n ^: 



rip) 



1 



7 



Dp 



2-7 



2- 7 

x1Z {2n) { Pl r) =5{p-r)p- {d -^ 

(d-2 + 4n)T(d - 2 + 2n) 
2 d - 1 Tr( d ~ 1 )/ 2 T((d- l)/2 + 2n) 



2k 



2n(2n + d-2) 



(3.6) 



where C^jf 1 was introduced in ( 1.14b ). 

Solving (3.€) one gets the behavior ofTZ^ 2 ^ in the dif- 



fusive and convective intervals explicitly: 

, (2n) _ (d-2 + An)T(d-2 



r, p, VTk <C r d 



2n) 



x exp 



Id/2-1 



2<*+l 7r (<2-l)/2 r ( 2n _ l + d/2) 
rp\ i \ l-d/2 



\rp 



8tK 

r,p, (Dt) 1 ^ »r d 1Z^ = 

(2 - 7)(d - 2 + An)T{d - 2 + 2n) 
7(d - l)2 d - 1 7r( d - 1 )/ 2 r(2n - 1 + d/2) 

(r-r + p7)( 2 - 7) \ /2(2- 7 )M^ 



(3.7a) 
(3.7b) 



x exp 



7 2 (d- l)W 



j 2 (d-l)tD 
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where 



1 /,, 8(d+l- 7 )n(2n + d-2) , . 
Vn = - ] J(d-^ + ^ ^ (3.8) 



and I v designate modified Bessel functions. Note that 
ft(o) = n { ^T{d/2)r 1 ' d /{2TT d / 2 ) where ft<» was intro- 
duced by ( |2.15U2.16D . It is pos sible to formulate for VS 2 ^ 
interpolation formulas of the (2.17) type. 

To establish the behavior of T n determined by 
we should analyze the integral (|3.3|) using (3.7). The 



main part of the integral is determined by the region 
r ~ ri3,r23. Thus the p -dependence is associated with 
the limit of (3.7) at p -C r. We see that the resolvents 



in this case tend to zero with decreasing p. We can also 
assert based on ([T|) that the higher is the number n 
of the angular harmonic the faster is the dec ay o f the 
corresponding resolvent. Since the procedure (|3 . 2|) pro- 
duces a convergent series over a small parameter then 
the above statements can be extended to the sum of the 
series. Thus we proved that the p-dependent part of T 
is paramctrically smaller than its p-independent part at 
p -C r\,r 2 what was the purpose of this subsection. 



B. Four-Point Correlator with two separated pairs 

At the rest of this section, we consider 1^1234 in the case 
of the special geometry with two separations between 
points being much smaller than all other ones, namely 
R > pi 2, where 



ri r 2 =pi, r 3 -r 4 = /9 2 , 
R= (r 3 + r 4 -r! -r 2 )/2. 



(3.9) 



Here, the mutual orientations of the vectors and the ratio 
pxj pi are supposed to be arbitrary. We denote by T (i?) 
the main contribution to 1^1234 which is equal to I1234 
at p\ = P2 = 0. The estimate for Tq can be extracted 
from the results that have been obtained in Section and 
justified in the preceding subsection. Indeed, the state- 
ments ( |2.54 ,2.57) are valid for r being the largest length 
R of the tetrahedron. One thus gets for To - . 

U < R < L T « (C1L 27 - C 2 R 2 ~<{L/R) A )PZ/D 2 , 

(3.10a) 

R^r d T « (C^ - C 3 RHL/R) A r 2 ;<- 4 )Pi/D 2 , 

(3.10b) 

where C\ , C 2 , C3 are dimcnsionlcss constants depending 
on d and 7. 

Generally, we do not assume here and below in this 
Section that A is small (or d is large) yet the condition 
A < 27 is assumed. At the limit d — > 00, the numerical 
coefficients C\, C2, C3 could be found from the results of 
Sect. II. In particular, we can get the following asymp- 
totic expressions: 



24 ( 2 -7) 2 ; , 2 , 

7 2 d 4 



(1 - 2/3) 



r 12. 



-2/3 



ri2 



Ai 



(3.11) 



at r 1 2«r 13 «r 2 3 ((0i - ~ ((^l ~ 

J13J \ria. 



4(2- 7 )>7 2 r7 3 



7 2 d 4 



(3.12) 



at 772, r 34 < r 13 w r 23 w r 24 ~ r 14 

« m - o 2 j 2 (e 3 - e 4 ) 2 ) - ((9, ~ e 2 f)((e z - e 4 f) 

4(2 - 7 ) 2 



7 2 (d- Ifd 2 



r 12 r 34 



L 

ri3 



(3.13) 



where dots designate the terms with a normal scaling 
(r 27 ). Note that if one expands ( 3.11 ) in small A-s then 
the remarkable cancelation of the terms linear in logs 
happens: (1 - 2/3)A 2 + 2/3Ai = 0. 

Now, we aim at finding leading corrections to Tq (R) at 
small values of pi/R and p 2 / R- For this we should solve 
fl2.6a ) at the conditions pi,p 2 "C R- In this limit the 
pri ncip al part of L is L Pl + C P2 where L p is determined 
by (3.1). It is the reason which caused us to formulate an 
iteration procedure where C Pl + C P2 is treated as a main 
contribution to C. Namely a solution of the equation 
( |2. 6a ) can be represented as 



ri234 = r (i?) + ^r„(R; Pl ,p 2 ) 



n=l 



-(4 1 +£ P2 )r„ = $„(R;p 1 ,p 2 ) 

$1 = /To + $12;34 + $13;24 + $14;23, 

$n = {C-L Pl -£pjr( n _i), 11>1. 



(3.14a) 

(3.14b) 
(3.14c) 
(3.14d) 



The procedure introduced by (3.14) can be considered 
as a series over the small parameters pi/R, p 2 /R- We 
impose two conditions on r„: to be of order of Tq at 
PI1P2 ~ R and to turn into zero at p\^ 2 = 0. 

Before consistent derivation, let us notice that the 
r^-dependence could be readily found from ( 3.1 4b| ) and 
the additional anomalous scaling of the subleading (p- 
dependent) terms naturally appears already at the first 
step as a consequence of the given scaling of Tq. Let 
us consider p <C ra, then the solution of (3.14b) for Ti 
will be p\ + p 2 multiplied by the i?-dependent terms 
$ and £Tq which both are proportional to H y (L/R) A 
- see (3.21) below. Assuming that no special cance- 
lations happen, one may conclude, in particular, that 
((0 2 (O)e(i?))) oc rJ~ 2 R^(L/R) A instead of R 2 ^- 2 (L/R) A 
that would be given by a naive counting of powers. 



Let us start now a regular examination of (3.14b). A 
solution of that equation has the following form 
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r„(R; Pl ,p 2 ) = dt (3.15) 
Jo 

x / / dr'dr"Tl(t; pi,r')TZ(t; p 2 ,r")$ n {K;r' ,r") , 



where TZ is the ke rnel of the resolvent of the operator 
C p introduced by (3.4,3.5 3.7). The main contribution 
to the integral (3.15) is determined by the region r' ~ R, 
r" ~ R. The higher is the number of the angular har- 
monic n m HQ the smaller is the respective con- 
tribution into r due to additional powers of the small 
parameter p/R in TZ^ 2n \ Therefore besides the first har- 
monic, we consider further only the second one to show 
how a non-trivial anomalous exponent appears due to an 
angular dependence. 



C. Zero modes 



The asymptotic formulas ( |3.7D for the resolvent enable 
one to extract a scaling behavior of the general expres- 
sion (3.15) for a solution of (3.141). We consider the 



case where the largest scale R of the tetrahedron is much 
smaller than the pumping scale L. That m eans that the 
only external scales for the integrations in ( 3.15 ) can be 
R or r d but not L. To be interested in the contribution 



to T n stemmed from large enough scales 



R and 



time t ~ tr (where one defines r r as r 2 / 'n at r <C r d but 
as r 1 / D at r > r d ), one can expand the resolvents in 
(3.15) in a series over the small parameters r r /t and ana- 
lyze the convergence of the temporal integral for different 
terms of the expansion. Of course, the integrals depend 
on the concrete form of <I>. However, in any case only a 
finite number of the first terms of the expansion are deter- 
mined by the time t ~ tr. The respective contributions 
into r are the zero modes of the operator C Pl + C P2 that 
grow with p and are symmetric under the permutation 
Pi *-* p 2 - All the rest terms (stemmed from r',r" <C R) 



are summed up into the term estimated as r p $, we will 
name this part of the solution the forced one 

The r-dependence of the forced part of the term r„ de - 
pends on a concrete form of the function $„ in (3.14b). 



Contrary, the form of the zero modes is universal: they 
depend on the source term only via coefficients. We 
thus start from finding out the zero modes that give 
the main contributions at small scales. To find the zero 
modes that contribute to the solution, we will follow an 
indirect (but, probably, the simplest) way. Since the zero 
modes are formed at the large times t ~ tr in the integral 
( pTT5| ) (which means that they appear due to matching 
at p r*> R). then we can e xtra ct their scaling behavior 
expanding the e xpre ssions (3.7) for the resolvents in the 



series over r p /t ( p.7|) and keeping the first terms of the 
expansion. Thus we can extract the exponents of the 
first zero modes and their angular dependence but not 
coefficients at different terms of the same order. To find 
the coefficients we can construct a combination of a given 
order with arbitrary coefficients and then find relations 
between the coefficients demanding that the combination 
is a zero mode of C Pl + C P2 , what can be established by 
direct applying the operator. 

Let us realize the scheme using the expressions ( [j.7| ). 
We see that the p-expansion of the expressions produces 
powers of p 7 and p s where as follows from (|3.i 



5=^,- d+ ,l( d - 1 ) 2 + 8{d+1 -" )d 



d-l 



(3.16) 



The angular structure of the zero modes also can be es- 
tablished if to take into account the explicit form of two 
first Jacobi polynomials: 



d + 1 



(dx 2 - 1) 



Then one can directly check that the zero modes of the 
first and of the second orders are as follows: 




P1P2 ~ 2pT7)(Pi 7 +> 9 2 7 )' P^> r d, 



(3.17a) 



{{pi~ 2 (d(Rp 1 ) 2 -R 2 p 2 )} + { Pl ^ p 2 })R- 2 , p>r d , 

{{d{Kp 1 ) 2 -R 2 p 2 } + {p 1 ^p 2 })r s d ~ 2 R- 2 1 p«r d ; ^' Ut>) 



~i"^> I£jrr-S(±\ J {(P2-M>? 2 (d(K Pl ) 2 - R 2 p 2 )}R- 2 + { Pl ^ p 2 }, p»r d , 
~ 2 ' 1 r>] ' W 2 -b 2 (P 1 ){d{^p 1 f-RV^}R- 2 + {px^P2}]rT S -\ P«r- d ; 1 ' 

Et R 2 h -8)(L\ A l ({^ 2 (d(R Pl ) 2 -i?V?)}{Pi-P2})i?- 4 + C i(piP 2 ) 5 - 2 (d(PiP 2 ) 2 -^), P»r d , 
2 ^ 2 \RJ \ [(KRpi) 2 - J RVf}{pi-P2})i2- 4 +c 2 (d(p lP2 ) 2 -p?pI)]rf- 4 , p«r d ; 

(3.17d) 
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where c\.i are dimensionless constants and 
((5 + 7 )(d + 2 7 - 



2d 



c?7 



(3.18) 



The upper indices of Z, introduced in ( p.l7| ), character- 
ize the type of a zero mode with respect to the angu- 
lar structure. The lower index denotes the original or- 
der of the resolvent's expansion over T p /t producing the 
concrete term. We keep the terms of the second order 



(3.17a, 3.17c, 3. 17d) besides the first order term (3.171) 
since they possess qualitatively different p-dependence: 
the term (3.17b) is an additive one while the terms 
(3.17a,3.17c,3.17d) contain cross-contributions. 

Note that higher-order angular harmonics will be char- 
acterized by the exponents 5 n which can be ext racted 
from the asymptotic behavior of the resolvents (3/7) char- 
acterized by ([T^)- Let us write the explicit expressions 
for the exponents: 



1 „ 8n(d+l- 1 )(2n + d-2) \ 
S n = ^[7-d+ ] J(d- 7 y + _ j 

(3.19) 



which at n = 1 gives ( |3.16 ). The exponents 5„ will figure 
in the expressions of the (3.17) type for higher harmonics. 
It is worthwhile to emphasize here that the exponents 6 n 
are nontrivial: to obtain them by £^" 1 £i-expansion that 
we used in Section^: one would get terms like p 7 ln(p/i?) 
that should be summed up into p s . 

We introduced in (3.17) the dimensional factors (in 
front of the braces) with which the zero modes appear 
at r. The factors in front of the braces are taken from 



the analysis of the integral (3.15); since the contribu- 
tions stem from the region p ~ R (where the resolvent is 
not precisely known) then the dimensionless coefficients 
(including c\ and ca) could not be found within our ap- 
proach with one exception that will be described below 



in Subsection 3 E , The dependence of the zero modes on 
rd can be established by matching the respective term 
from inertial and diffusive interval at and using the 
rrf-independence of T in the convective interval. 



D. Forced Terms 



Let us begin the discussion of the forced terms with 
r{ that appears at the first step of the iteration proce- 



dure (3.14c). The r.h.s. of (3.14c) contains the terms 
of different order in p/R, different angular functions and 
different types of the dependence on pi and p 2 (addi- 
tive and multiplicative). We shall analyze the respective 
contributions order by order in p/R separately for the 
different angular functions. The terms with the cross- 
dependence (like P1P2) are °f special interest since they 
contribute into (ee) correlation function. 



The leading in p/R contribution in r.h.s. of ( |3.14c ) is 
p-independent (and consequently angular independent). 
The contribution originates 
a) from CTq: 



CT ^(£ r -kA r )T (R); 
b) from <J>i 3 ;24 + $14;23: 



(3.20a) 



$13.24 + <J>14;23 - *R = 2 ^^(^) 2 - ( 3 - 20b ) 



The expression for the respective forced solution is 
(2-7) 



l.add 



(d-l)dD 
rdr 



<f>n + (C R - kAr)To(R) 

P2 



+ r 2 ~ 7 







2r 2 " 7 



+ r 2 -r 



(3.21) 



The subsequent additive terms are produced by taking 
into account the next p 2 -terms at the r.h.s. of (3.14c). 
They produce the additive forced contributions to Ti 
with the scaling oc p 2+7 in the convective interval. 

The leading cross-term originates from $12.34 , the 
main term of its expansion is 



g P?P?^d/(pi)d/( P2 ) 
P1P2 dpi dp 2 



12;34 



(3.22) 



Apart from (3.22), we have to take into account the cross- 
terms originating from other terms in r.h.s. of (3.14c). 
The terms have a character of a regular expansion in p/R 
and are consequentl y pro portional to p\p\. Such terms 
are much less than (3.22) in the convective interval and 
are of the same order in the diffusive interval. Since we 
are interested not in the factors but only in the scaling 
behavior we can restrict ourselves with the analysis of 
( |3.22 ) . The term is proportional to p 27 in the convective 
interval and to p 4 in the diffusive one, that gives for the 
principal behavior of the corresponding forced term p 37 
and p 6 respec tively . 

The source ( 3.22 ) produces not only cross-terms in L of 
scalar and tensor structures but also additive and gener- 
ally mixed ones possessing all the variety of angular and 
Pi P2 symmetries used in the classification of the zero 
modes ( 3.17 ). All those terms scale with R as i?~ 7 both 
in the convective and diffusive intervals by p. To match 
those forced terms from the convective and diffusive in- 
tervals one has to introduce the additional zero modes in 
the diffusive interval by p 



cr(0,0) 
-°2,d 


^Z^(R/L) A 


[u/Rf, 


(3.23a) 


7(0.2) 


^Z[°' 2) (R/L) A 


{r d /Rf~<-\ 


(3.23b) 


7(0.2) 
-°2,d 


^^(R/Lf 


[rd/R) 2 ~<-\ 


(3.23c) 


7(2,2) 


^Z^(R/Lf 


{rd/Rf- 26 . 


(3.23d) 
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Let us now consider the forced term arising at the 
second step of the iteration procedure. Analyzing the 
structure of the r.h.s. of (3.14d) one concludes that it 



decreases downscales with p not slower than (r{+zcro 
modes). Consequently, the forced term r| decreases 
downscales with p not slower than T p (r{+zero modes) 
that is faster than r{. Thus, one can drop in com- 
parison with r{ and all the set of zero modes. The same 
arguments allows one to drop the higher forced terms r£, 
n > 2 too. 



E. Comparison between the zero modes and forced 
terms 



in the convective and diffusive intervals respectively. The 
zero modes ( 3.17b| , 3.17c ) and ( 3.17c] ) scale as p 5 , p s+1 



and p 2S in the convective interval, and as p , p and p 
in the diffusive one respectively. One concludes: 

a) in the convective interval, the zero mode {zf^ 

(2 2) 

or ' ) is dominant if 37 > S (respectively 2j > 5 or 
7 > 2i5/3), otherwise at the smallest p the forced terms 
prevail; 

b) in the diffusiv e inter val, at 37 > S (2 -f > 6 or 7 > 26/3 ) 
the zero mode ( 3.17b ) [respectively ( 3.17c ) or ( 3.1 7d )] 
becom e dom inant at small enough p; otherw ise the zero 
mode ( 3.23b ) [respectively ( 3.23c ) or ( 3.23d )], matching 
with the forced solution originating from (3.22), prevails. 



To find the main contributions into T one should com- 
pare the zero modes and forced terms of the same angular 
structure. The leading term in the expansion in p/ R is 
the isotropic additive forced term (3.21) inside of both 
the convective and diffusive intervals. 

The zero modes are L-dependent but the non-additive 
forced terms are not. Thus keeping the parameter L/ R to 
be large enough one forces the zero modes to dominate. 
A possibility for the forced terms to prevail appears at the 
smallest p-separations, when the dominance of the zero 
modes could be compensated by different p-dependencies 
of the zero modes and the forced terms. Present Subsec- 
tion is devoted to such a comparison. 

We co mpare ( 3.17a ) with the forced terms generated 
by ( 3.22 ). Both in the convective and diffusion interval 
we get the d ominance of the angular independent zero 
mode ( 3.17a) co ntaining the cross-contribution. To be 
sure that ( 3.17a ) does give the main contribution into 
{{^1^2)) (and thus provides for its anomalous scaling), we 
have to check whether or not the respective dimensionless 
coefficient in front of that mode can turn into zero. We 
consider that mode in a special case R L p when 
the coefficient can be found unambiguously contrary to 
other cases where the matching at p ~ R is necessary. 
If R 3> L ^ p then one have to consider the integral 
(3.15) with (3.22) and the respective term from (3.7) of 
the expansion of - see (2.16b 2.17). One can directly 
check that the main contribution into the integral stems 
from the scales L and time tl so that no matching at R 
is necessary. We do not write here the bulky expression 
for the integral, it turns into zero only at a single value 
"f/d w 0.086 found numerically. That numerical value 
depends on the shape of x- At the rest values of j/d the 
integral is nonze ro and provides the contribution of the 
structure ( |3.17a| ). 

Next we consider the second angular harmonic pre- 
sented in (3.17b ,3.23b), (3. 17c. 3. 23c) and the product of 
the second angular harmonics in ( 3. 17d, 3.23d ) . The scal- 
ing of the respective harmonics of the forced terms coin- 
cides with that of the isotropic contribution: p 37 and p 6 



4. POSSIBLE GENERALIZATIONS 

Let us remind the reader that the above results on 
r^-dependence are general while those on L-dependence 
are formally obtained only when the respective anoma- 
lous exponent A is much less than the normal exponent 
27. The whole approach of Sect. 2 was developed for the 
fourth-order correlation function at d > 2. The general- 
ization for the n-th correlation function at d > n — 2 is 
straightforward by means of the same representation of 
C in terms of . Such a generalization leads to the new 
predictions for n > 4: the scaling exponent 27— A should 
appear at all high-order functions. The direct check 
shows that the respective zero modes appear in the par- 
tially reducible contributions only. For example, the six- 
point object F 123456 = (0(ri)0(r 2 )e(r 3 )e(r4)0(r 5 )0(r 6 )) 
should necessarily contain terms like /(?"i2)r3 4 56 i.e. 
r 27 ~ A L 7+A . Such contributions should cancel out in 
structure functions like (((#1 — #2)™)). To find the scal- 
ing exponents of the n-th structure functions in the limit 
d — > 00, one should iterate the zero mode with the scaling 
exponent 717, which will be published elsewhere. 

Quite a different picture may appear for d < n — 2. 
In this case, n(n — l)/2 distances between points are not 
longer the independent variables. Additional constraints 
should be imposed on that set of (otherwise very con- 
venient) variables, which may lead to the possibility of 
additional zero modes. 

As a final remark, let us emphasize that the above 
results are valid for a Gaussian pumping only and dis- 
cuss what might be the consequences of the pump- 
ing non-Gaussianity. Let us add to the fourth-order 
pumping correlation function ( |1.4| ) the irreducible part 
X4(ri,r 2 ,r3,r 4 )5(ii - t 2 )S(t 1 - t 3 )(5(ii - i 4 ). One 
may model the function % 4 by the step function: 
X 4 l (ri, r2, r3, r 4 ) = P 4 if all 77, < L and zero if any 
rij > L. The production rate of 8 A is 3P2/(0) + P 4 (let us 
remind that the advection preserves the integral of any 
power of the scalar field). The ratio r* = P4/P2 hav- 
ing the dimensionality of time is a correlation time of the 
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pumping, which is assumed to be the smallest time in the 
problem. As a result, the terms generated by Xi wm be 
small yet we keep them to show what qualitatively new 
may appear in the case of a non-Gaussian pumping with 
a finite correlation time. The equation ( |2.6a ) acquires 
additional term: 

— ^ri234 = Xi + ^12:34 + < f ) 13;24 + $14;23- 

The presence of the function Xi changes scaling in the 
convective interval, it causes the terms with anomalous 
scaling T*r^- to appear in T which corresponds to the 
scaling ( n = ( 2 discussed in Jl^,§|. Also, the terms with 
r® 2 appear in ((ei0|)), while no additional terms appear 
in ((£162)). Note that within our approach the r„-related 
terms are small corrections even at the diffusion scale 
despite the fact that they decrease slower than the main 
terms as distances decrease. Most probably, finite cor- 
relation time of the pumping will lead to a substantial 
contributions with an another anomalous scaling yet this 
is beyond the scope of the present analysis. The develop- 
ment of the theory for the finite correlation times of the 
velocity and pumping is necessary for a meaningful com- 
parison between theory and experiment. This will be the 
subject of future publications. On the other hand, the 
more detailed experimental results are desirable, in par- 
ticular, direct measurements of the correlation function 
of the dissipation field. 



5. CONCLUSION 

We have shown that the fourth-order correlation func- 
tion F1234 has the scaling exponent 27 — A in the con- 
vective interval at d > 2. The anomalous exponent A is 
analytically found at d — > 00. Considering the behavior 
of -F1234 in the convective interval where some separa- 
tions p tend to zero we established that the p-dependent 
contribution to F1234 also tends to zero. Moreover, noth- 
ing terrible happens when p passes through the diffusion 
scale r c i and we find that 



9 2 ) 4 ) oc |n 



T2 



27-A^A 



(5.1) 



for r*i2 from the convective interval. 

We analyzed also the case where two separations are 
much smaller than the distance R between the pairs 
and found the contribution, mixing both small distances, 
which determines the correlator 

((01 - 2 ) 2 (9 3 - 8 A ) 2 ) - ((9, - 9 2 ) 2 ) ((9 3 - 9 4 ) 2 ) , (5.2) 

at rn, ^34 <C R- We established that the scaling be- 
havior of that correlation function in the convective in- 
terval is characterized by the additional anomalous ex- 



ponent 5 defined by (3.16). The corresponding contri- 
butions to the irreducible part 1^1234 of the correlator 



-F1234 are presented in (3.17). The isotropic contribution 
( |3.17a ) behaves oc rJ 2 rJ 4 (L / R) A , while the anisotropic 



one ( |3 17d|) b ehaves 

^ r i2 r 34^ 2S \L/R) A . The cross 



term ( |3.17c| ) behaves oc (rf 2 r^ 4 + r d 34 rJ 2 )R^-° (L/ R) 1 
The analysis of the continuation of the correlation func- 



tion (5.2) for 7*12, passing to the diffusive interval 
gives the following scaling laws 

«e ie3 » = k 2 «(V0i) 2 (V0 3 ) 2 » « r d (L/r 13 ) A , (5.3) 
(((V^iV^i - (TX/sCWi) 2 ) 

x(V„0 3 V„0 3 - d-H^{V9 3 ) 2 ))) oc r 2 d s - 4 r 2 r 25 (L/r 13 ) A , 

(5.4) 

((eilV^V^s - d" 1 ^ (V0 3 ) 2 ]» oc r s - 2 rJ~ S (L/r 13 ) A , 

(5.5) 

where the separation n 3 lies in the convective interval. 
The double angular brackets designate irreducible cor- 
relation functions. The analogous analysis can be per- 
formed for the correlation function 



((ei-9 2 ) 2 el)-((e 1 -e 2 f)(el) 



(5.6) 



at 7*12 <C ri 3 . The contribution (3.21) gives the law 
oc v1 3 rj 2 for the isotropic part of the correlation func- 
tion (|5.6|) what leads to 



((<*%)) <xrfrl 3 {L/r 13 ) 



(5.7) 



The anisotropic part of the correlation function (5.6) is 
determined by the contribution (3.17b) what is propor- 



tional to r 2 g S rf 2 , it leads to 

{{el[V a 8 2 Vp8 2 - d- 1 S al3 (V8 2 ) 2 ])) oc r s d - 2 rlr 5 {L/r l3 ) A . 

(5.8) 

The ap pearance of r^-dependence in the correlation 
functions (5.3-5.8) may correspond to the presence of 
the ultraviolet divergences found in jl6|,|ll| in the dia- 
grams for powers of gradients both of the velocity and 
of the passive scalar. That means that, contrary to L- 
related scaling, r^-related anomalous scaling of the scalar 
derivatives could be caught perturbatively. Moreover, 
in the asymptotic limits of the normal overall scaling 
of Ti234 (at d — > 00 or 7 — > 2 ~) it is tempting to de- 
scribe the structure of (5.1,5.3,5.4.5~5|,5.7 5~8|) u sing the 
language of the so-called operator algebra |)26|-p8| (de- 
veloped in the context of second-order phase transition 
theory) the validity of which for turbulence was argued 
in jui]]. Namely, at d — > 00 or 7 — * 2~ in accordance 
with those expressions, the passive scalar 8 has the di- 
mensionality —7/2, (V6 1 ) 2 has the dimensionality and 
(V a 9V j39~d~ l 8 a j3(V9) 2 ) has the dimensionality <5— 7. In 
the general case of arbitrary d and 7 when the L-related 
anomalous scaling does exist, the underlying algebraic 
structure (if it exists at all) is unclear at the moment. 
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All said above concerns the terms associated with 
zero modes of the operator C in the convective inter- 
val. There exist also a forced solution in the convective 
interval originating from ( 3.22| ). The tails in the diffu- 
sion intervals created by this forced solution can be in- 
terpreted as zero modes (in the diffusive interval only) 
which are determined by the estimates (3.23). Those 
zero modes produce the contributions to the correlation 
functions (5.4 jp5 |5~8| ) which are oc r^(L/r1 3 ) and pre- 
vail if S > 37/2, S > 7/2 and S > 3j correspondingly. 
Those contributions are not governed by the operator al- 
gebra even in the asymptotic cases of d — > 00 or 7 ^ 2~. 
It is natural from the diagrammatic point of view since 
those contributions correspond to "one-bridge" diagrams 
which does not contain series producing anomalous scal- 



ing. Thus the factor 



13 



is determined simply by the 



"bridge" factor VV/C oc R~~> . Of course analogous con- 
tributions exist in all above correlation functions but only 
in (5.4,5.5,5.8) they can prevail at some 5. 



To conclude, 1/d-expansion tells us that the anoma- 
lous scaling is present already in the oversimplified model 
of the delta-correlated velocity field. If one formally 



uses our formula (3.1) for the anomalous exponent at 
Richardson-Kolmogorov case d = 3,7 = 2/3 it gives 
A = 16/9 which substantially larger than the experi- 
ments and numerics give (A ~ 0.4 0.5), the difference 
could be accounted by inaccuracy of both an asymptotic 
formula at finite d and the model with delta-correlated 
velocity. Note that our results are valid not only for 
steady state but also for decaying turbulence of a pas- 
sive scalar with a compact spectrum as an initial condi- 
tion. In that case, the scalar cascade is accelerated as it 
goes towards small scales (with typical time t oc r^ 2-7 " 2 ) 
so that the small-scale part of the scalar distribution is 
quasi-steady. 
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